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Abstract 

We derive an energy-independent and non-local potential for the baryon-baryon interaction from 
the Nambu-Bethe-Salpeter amplitude on the lattice. The method is applied to the nucleon-nucleon 
interactions with the (2-(-l)-flavor lattice QCD simulations. The central and tensor potentials are 
obtained as leading order terms of the velocity expansion of the non-local potentialF The central 
potential has a repulsive core surrounded by an attractive well, while the tensor potential has inter- 
' mediate and long range attractions. Validity of the velocity expansion is tested by the nucleon-nucleon 

interaction with finite relative momentum. Interactions among octet-baryons in the flavor SU(3) limit 
are also studied in full QCD simulations to unravel the origin of the repulsive core and the possible 
existence of _ff-dibaryon. 
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1 Introduction 



Understanding of the nuclear force from quantum chromodynamics (QCD) is one of the most challenging 
problems in nuclear and particle physics. Experimentally, a large number of proton-proton and neutron- 
^ ' proton scattering data as well as deuteron properties have been accumulated and summarized e.g. in 

P5 , the Nijmegen database ^1 . Below the pion production threshold, the notion of the NN potential (either 

in the coordinate space or in the momentum space) is useful in the sense that it can be used not only 
to describe the two-body system but also to study the nuclear many-body problems through ab-initio 
^ ' calculations [2]. 

CO i The phenomenological NN potentials in the coordinate space are known to reflect some characteristic 

' features of the NN interaction ;3]: 

(i) The long range part of the nuclear force (the relative distance r > 2 fm) is dominated by the one-pion 
exchange [4. Because of the pion's Nambu-Goldstone character, it couples to the spin-isospin density of 
the nucleon and hence leads to the tensor force. 

(ii) The medium range part (1 fm < r < 2 fm) receives significant contributions from the exchange of 
multi-pions. In particular, the spin-isospin independent attraction of about 50 - 100 MeV in this region 
plays an essential role to nuclear binding. 

(iii) The short range part (r < 1 fm) is best described by a strong repulsive core [5l l6l . Such a short 
' range repulsion is relevant for the maximum mass of neutron stars. 

^ I (iv) There is also a strong attractive spin-orbit force in the isospin 1 channel at medium and short 

distances. This leads to the neutron pairing in neutron matter and hence the neutron superfluidity 
inside neutron stars JT^. 

Several high precision NN forces are now available to fit neutron-proton and proton-proton scattering 
data (about 4500 data points) with x^/dof ^ 1. However, they have typically 20-40 fitting parameters: 
e.g. CD Bonn potential, AV18 potential and N'^LO chiral effective field theory have 38, 40, and 24 param- 
eters, respectively 8 . If one tries to extend these to hyperon-nucleon and hyperon-hyperon interactions, 
the task becomes extremely tough since the number of parameters increase and the scattering data are 
scarce. In this situation, it is highly desirable to study the general baryon-baryon interactions from the 
first principle lattice QCD simulations, since all the hadronic interactions in QCD are controlled only 
by the QCD scale parameter (Aqcd) and the quark masses (to„, m^, m^) whose values are pretty well 
determined [9j. 

A theoretical framework to study the hadron-hadron interaction using lattice QCD was first proposed 
by Liischer [TU] and was applied to the lattice simulations for the NN interaction in [IT: For two hadrons 
in a finite box with a size L x L x L in the periodic boundary condition, an exact relation between the 
energy spectra in the box and the elastic scattering phase shift at these energies was derived. If the 
range of the hadronic interaction R is sufficiently smaller than the size of the box R < L/2, the behavior 
of the equal-time Nambu-Bcthe-Salpeter (NBS) amplitude il'{r) in the interval R < \r\ < L/2 under 
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the periodic boundary condition has sufficient information to relate the phase shift and the two-particle 
spectrum. This Liischer's method bypasses the difficulty to treat the real-time scattering process on the 
Euclidean lattice. 

Recently, an alternative approach to the hadron interactions in lattice QCD was proposed p!^[T51[H] . 
The starting point is the same equal-time NBS amplitude tp{r): Instead of looking at the amplitude 
outside the range of the interaction, the internal region |r| < i? is considered and an energy-independent 
non-local potential U{r,r') is defined from ip{r). Since U{r,r') in QCD is a localized function in space 
due to confinement of quarks and gluons, it receives finite volume effect only weakly. Therefore, once 
?7, although it is not a direct physical observable, is determined on the lattice, one may simply use the 
Schrodinger equation in the infinite space to calculate observables such as the scattering phase shifts, 
bound state spectra etc. Moreover, the potential would be a smooth function of the quark masses: This 
is in sharp contrast to the scattering length which shows a singular behavior around the quark mass 
corresponding to the formation of the two-body bound states such as the deuteron. Similar situation is 
well-known in the BEC-BCS crossover of cold fermionic atoms jl5j . 

In this article, we will show some recent results of the nuclear force on the lattice (or the lattice 
nuclear force in short) after a brief introduction to our basic formulation. 

2 Deriving the NN potential on the lattice 

2.1 NBS wave function on the lattice 

In field theory, the best analogue of the two-particle wave function is the equal-time Nambu-Bethe- 
Salpeter (NBS) amplitude or the "NBS wave function": Let us consider an exact six-quark state \£) 
which has total energy £, total three-momentum zero and total electric charge +e in a finite box. Then 
we define the NBS wave function by 

V'(r) = {0\np{x + r,t = 0)p^{x,t = 0)\£). (1) 

The local composite operators for the proton and the neutron are denoted by Pa{x,t) and np{y^t) with 
spinor indices a and /3. One should keep in mind that \£) is not a simple superposition of a product state 
p) (g) |n), since there are complicated exchanges of quarks and gluons between the two composite particles. 
The NBS wave function ^{r) can be regarded as a probability amplitude in \£) to find "neutron-like" 
three-quarks located at point x + r and "proton-like" three-quarks located at point x. 

The spatial extent of the NN interaction in QCD is short ranged and is exponentially suppressed 
beyond the distance R ^ 2 fm. Therefore, the spatial part of the NBS wave function in the "outer region" 
satisfies the Helmholtz equation, 

(V^ + fc2)V.(r) = (|r|>7?), (2) 

up to an exponentially small correction in L. Here the "asymptotic momentum" k is determined by the 
asymptotic behavior of the wave function in the outer region. 

An important property of the NBS wave function ip{r) is that its asymptotic behavior at large \r\ in 
the infinite volume limit reproduces the correct phase shift obtained from the iS-matrix of the elastic NN 
scattering. This can be shown explicitly by using the Nishijima-Zimmermann-Haag(NHZ)'s reduction 
formula [16J for the products of local composite operators. (See Appendix A of [14 for the proof.) To 
define the NBS wave function on the lattice, we start with the four-point function 

g{r,t^to) = {0\np{x + r,t)pa{x,t)S{to)\0) ^ ^/j{r) e-^o(*-*o) » to), (3) 

where £q is the lowest energy state created by the source operator 5(io)- 

2.2 Non-local potential and the velocity expansion 

To define the NN potential from the NBS wave function, let us introduce the following local function: 

KE{r) = ^(V^ + fc2)^Jr) = {E- Ho)i>,{r). (4) 

In the second equality, we introduce an "effective center of mass energy", E = k'^/{2fi), and the free 
Hamiltonian Hq — — V^/(2/i), with /i = mN/2 being the reduced mass of the two nucleons. They are 
introduced only to make a formal resemblance with the Schrodinger type equation and have nothing 
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to do with non-relativistic approximation. Hereafter, we put the suffix E to the NBS wave function to 
emphasize its i?-dependence. Since the "plane-wave" part of the NBS wave function in the outer region 
(r > R) is projected away by the operator E — Hq, the function KE{r) is non- vanishing only in the inner 
region (r < R). Note also that the Fourier transform of Ksir) is essentially the half off-shell T-matrix. 
We can rewrite Eq.Q in two equivalent ways: 

iE-Ho)iPAr) = UE{r)i;^{r)= f U{r,r')ij^{r')dr' . (5) 



The first equality is just a definition of the energy- dependent local potential, Ue{t') — K E{r) / j^{r) . On 
the other hand, the energy-independent non-local potential, U{r,r'), is defined from Ue{t') through a 
self-consistent equation, 

U{r,r') = {r\U\r') = V / UEir){r\e'^""+^~^^*\r'). (6) 



E 

Carrying out the t integration formally, one may also write Eq.® as C/ = UeS{E — Hq — U). In 
these formulas, stands for the summation (integration) over the discrete (continuum) energies. In 
particular, E is always discrete on the lattice with a finite volume. Also, E has an upper limit Ec at which 
inelastic scattering starts to take place. Eliminating the i?-dependence of the potential through Eq.([B]) has 
been discussed in a transparent manner by Krolikowski and Rzewuski [17] long time ago: their motivation 
was to prove the equivalence between the multiple-time Nambu-Bethe-Salpeter type equation with an 
-E-dependent kernel and the equal-time Schrodinger type equation with an £'-independent potential. 
Essentially the same method was rediscovered and discussed in [121 HSl E] in the context of the NBS 
wave function on the lattice. 

If we further focus on the low-energy scattering with E sufficiently smaller than the intrinsic scale of 
the system or the scale of the non-locality of the potential in Eq.®, the velocity expansion of U{r, r') in 
terms of its non-locality is useful |18| : For example, the potential with hermiticity, rotational invariance, 
parity symmetry, and time-reversal invariance may be expanded as |19] 

[/(r,r') = V{r,v)5{r -r'), (7) 
V{r,v) = Vcir) + ^7(^)^12 + VLs{r)L ■ S + 0{v^) + ■■■, (8) 

LO NLO N^LO 

where v = p/^ and L = rxp with p — and 5*12 = 3{ai-r){(T2-r)/r^ — (Ti-(T2- Each coefficient of the 
expansion is a local potential and can be determined successively by measuring the NBS wave functions 
for several different energies. The central potential Vc and the tensor potential Vt are classified as the 
leading order (LO) potentials since they are of 0(v°). The next-to-leading (NLO) potential of 0{v) is 
the spin-orbit potential VLs(r). The LO and NLO potentials are phenomenologically known to be the 
dominant interactions at low energies. 

An advantage of defining the potential from the NBS wave functions in the "inner region" is that the 
effect of the lattice boundary is exponentially suppressed for finite range interactions: Then one can first 
make appropriate extrapolation of U{r, r') or V(r, v) to L ^ 00, and then solve the Schrodinger equation 
using the extrapolated potential to calculate the observables such as the phase shifts and binding energies 
in the infinite volume. This is in contrast to the Liischer's approach |10| in which the wave functions in 
the "outer region" suffering from the boundary conditions is ingeniously utilized to probe the scattering 
observables. Apparently, the two approaches are the opposite sides of a same coin. 



2.3 Interpolating operator and the potential 

In Eq. ([1]), simplest interpolating operators for the neutron and the proton written in terms of the 
up-quark u{x) and the down-quark d{x) would be 

ni3{x) = Eabc {Ua{x)C'-f^db{x)) dcp{x), Pa{x) = £abc {Ua{x)C-f5db{x)) Uca{x), (9) 

where x — {x,t) and the color indices are denoted by a, b and c. The charge conjugation matrix in the 
spinor space is denoted by C. The local operators given above are most convenient for relating the NBS 
wave function to the four-point Green's function and the scattering observables at L oo through the 
NZH reduction formula. 

In principle, one may choose any composite operators with the same quantum numbers as the nucleon 
to define the NBS wave function. Different interpolating operators lead to different NBS wave functions 
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Figure 1: (a) LO potentials in (2+l)-flavor QCD for m^=701 MeV [H]. (b) LO potentials in quenched 
QCD for m^=731 MeV [M]. 
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Figure 2: Quark mass dependence of the LO potentials in (2+l)-flavor QCD. (a) The central potential in 
the spin-singlet channel, (b) the central potential in the spin-triplet channel, and (c) the tensor potential 
in the spin-triplet channel [22]. 



and different NN potentials. However, they lead to the same physical observables by construction. 
Analogous situation can be seen in quantum mechanics where the unitary transformations modify both 
the wave function and the potential in such a way that observables are unchanged. Even more direct 
analogy is in field theory for point-like particles: Field re-definitions modify the vertices and propagators 
in the Feynmann rule, while the on-shell S'-matrix is not affected by such changes. 



2.4 Central and tensor forces 

In the LO of the velocity expansion in Eq. ([5]) , we have the central potential Vc (r) and the tensor potential 
Vrir), so that the Shrodinger equation reads 

{E - Ho)i^,{r) = {Vc{r) + VT{r)S,2)i^,{r). (10) 

The central potential acts separately on the orbital S-state and the D-state, while the tensor potential 
provides a coupling between these two. Therefore, a coupled-channel Schrodinger equation is obtained 
from Eq. (fTO|) by using the projection operators P and Q to the S-state and D-state, respectively. Even- 
tually we calculate Vc and Vr from the following formula where the quantities in the right hand side are 
all known on the lattice: 

( Vc\_f Vi^, VSi2i^, ( E-Ho V ^ nn 

{Vt J \ Q^, QS,2i', J \ E - Ho J \ Qi;, J ■ 



3 Numerical results in quenched and full QCD simulations 
3.1 LO potentials 

To demonstrate whether the formalism discussed in the previous section indeed works, we first carried 
out a quenched QCD simulations with the standard plaquette gauge action and the standard Wilson 
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Figure 3: (Left) A comparison of the central potentials obtained at different energies (E ^ MeV vs. 
i? ~ 45 MeV). (Right) A comparison of the spin-singlet central potential at i? MeV with different 
orbital angular momenta (L — vs. L — 2). Figures are taken from |23) . 

quark action on a 32'^ x 48 lattice [121 [131 [14] . The lattice spacing is a = 0.137 fm which corresponds to 
the spatial size L = 4.4 fm. The light quark masses are chosen so that we have =731, 529 and 380 
MeV and =1558, 1334 and 1197 MeV, respectively. Periodic or anti-periodic boundary conditions 
are imposed on the quark field along the spatial direction. 

As for full QCD with the vacuum polarization of quarks included, we use the PACS-CS gauge con- 
figurations in (2-|-l)-flavor QCD generated by the Iwasaki gauge action and the 0(a)-improved Wilson 
quark (clover) action on a 32'^ x 64 lattice [21]. The lattice spacing is a = 0.091 fm which corresponds to 
the spatial size L = 2.9 fm. The light quark masses are chosen so that we have TO,r =701, 570 and 411 
MeV and toa? =1583, 1412 and 1215 MeV, respectively. Also, Periodic boundary condition is imposed 
on the quark field along the spatial direction. 

Shown in Fig[TJa) are the LO potentials {Vc for ^5*0 and '^Si channels and Vr determined from ^Si- 
'^Di channel) in (2-|-l)-flavor QCD for to^=701 MeV. Even with such a large quark mass, there is a clear 
evidence of the repulsive core surrounded by attractive well for for central potential and an evidence of 
a mild tensor force [22]. They have qualitative similarity with phenomenological potentials. We show 
in Figlljb) the LO potentials in quenched QCD with m7r=731 MeV for comparison [2]. Although the 
qualitative structure of the potentials are the same, the magnitude of the repulsive core and the tensor 
force are relatively weak in quenched QCD. 

Shown in Figl2ja,b,c) are the quark mass dependence of the LO potentials in (2-|-l)-flavor QCD [22] . 
As the quark mass decreases, the repulsive core in (a,b) and the tensor force in (c) become stronger and 
the attractive well in (a,b) becomes wider. We have fitted these potentials and have calculated the NN 
scattering phase shift by solving the Schrodinger equation. We found that deuteron bound state does not 
appear for these quark masses, so that further reduction of the quark mass would be necessary to obtain 
the realistic lattice potentials. 

3.2 Convergence of the velocity expansion 

So far, the potentials are derived with the periodic boundary condition in the spatial direction for the 
quark fields. This leads to the the "effective center of mass energy" E = k'^ /{2fi) almost zero. To study 
the convergence of the velocity expansion of the non-local potential in Eq.®, we compare the local ^5*0 
potential (in quenched QCD with = 529 MeV) obtained at i5 ~ MeV under the periodic boundary 
condition and that obtained at ii^ ~ 45 MeV under the anti-periodic boundary condition [23 . Good 
agreement between the two as shown in the left panel of FigO indicates that a 0{v'^) term in the N^LO 
level is rather small in this energy interval. Shown in the right panel of FigO is a test for a different 
0{v^) term in the N^LO level [33]: In this case, local potentials determined at the same energy (i? ~ 
MeV) with different orbital angular momenta (L = 0, 2) in the spin-singlet channel are compared. Again, 
within statistical errors, the effect of the N^LO term is likely to be small. 
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Figure 4: NBS wave functions at = 835 MeV, normalized to 1/2 for the singlet channel and to 1 for 
other channels at the maximum distance |25|. 
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Figure 5: The three independent BB potentials in the channel in the flavor SU(3) limit, extracted 
from the lattice QCD simulation at TOtt — 1014 MeV (red bars) and m^r = 835 MeV (green crosses). 



4 Hyperon interactions 

To unravel the origin of the repulsive core in the NN interaction, let us consider the S'-wave interaction 
between octet baryons in the flavor SU(3) limit. In this case, two baryon states with a given angular 
momentum are labeled by the irreducible flavor multiplets as 

8®8 = 27©8^©J. ® 10* ® 10 ® 8a . (12) 
symmetric anti-symmetric 

Here "symmetric" and "anti-symmetric" stand for the symmetry under the flavor exchange of two baryons. 
For the system in the orbital S-wave, the Pauli principle between two baryons imposes 27, 8s and 1 to 
be spin singlet (^5o) while 10*, 10 and 8a to be spin triplet {'^Si). Since there are no mixings among 
different multiplets in the SU(3) limit, one can define the corresponding potentials as 

'So : V^^'Hr), T/(«=)(r), V^'\r), (13) 
35i : V^^°'\r), F(i°)(r), V'^^^\r) . (14) 

Potentials among octet baryons, both the diagonal part {B1B2 B1B2) and the off-diagonal part 
{B1B2 — > B3B4), are obtained by suitable combinations of F^"-'(r) with a = 27, 8^, 1, 10*, 10, 8q. 

In this SU(3) study, we employ the gauge configurations on a 16"^ x 32 lattice generated by CP- 
PACS and JLQCD Collaborations with the renormalization group improved Iwasaki gauge action and 
the non-perturbatively 0(a) improved Wilson quark action. The lattice spacing and the lattice volume 
are a = 0.121(2) fm and L = 1.93(3) fm, respectively. These configurations are provided by Japan Lattice 
Data Grid (JLDG) and International Lattice Data Grid (ILDG) [24] . 

Fig. m shows the NBS wave functions as a function of the relative distance between two baryons at 
rriTr = 835 MeV [25]. To draw all data in a same scale, they are normalized to 1/2 for the singlet channel 
and to 1 for other channels at the maximum distance. The wave functions in Fig. |4]show characteristic 
fiavor dependence: In particular, a strong suppression at short distance appears in the 8s channel, while 
a strong enhancement appears in the 1 channel. Similar results are obtained for 111^^ = 1014 McV too. 

Fig. [5] shows the resulting three independent BB potentials in the 'Sq channel in the flavor basis 
obtained from the NBS wave functions. Red bars (green crosses) data correspond to the pion mass 1014 
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MeV (835 MeV): Although there is a tendency that the magnitude (range) of the potentials becomes larger 
at short distances (longer at large distances) for lighter quark mass, the differences are not substantial 
for the present heavy quark masses. Left panels of Fig. [5] show which corresponds to NN ^5*0 

potential. It has a repulsive core at short distance and an attractive pocket as we have shown already 
in quenched and (2+l)-flavor simulations. The middle panel of Fig. [5] corresponding to has a very 

strong repulsive core among all channels. In contrast, the right panel of Fig. [5] corresponding to T/^-*-' 
shows attraction for all distances, which is relevant to the cereblated _ff-dibaryon [26] . 

These features are consistent with what has been observed in phenomenological quark models [27] . In 
particular, the potential in the 8s channel in quark models becomes strongly repulsive at short distance 
since the six quarks cannot occupy the same orbital state due to quark Pauli blocking. On the other hand, 
the potential in the 1 channel does not suffer from the quark Pauli blocking and can become attractive 
due to short range gluon exchange. Such an agreement between the lattice data and the phenomenological 
models indicates that the quark Pauli blocking plays an essential role for the repulsive core in BB systems 
as suggested long time ago in [28]. One can also confirm the idea of the Pauli blocking by considering 
the meson-baryon interaction such as charmonium-nucleon potential [29) and kaon-nucleon potential [30] 
within the present lattice approach. Generalization of the baryon-baryon interaction to the case with 
exphcit SU(3) breaking is also under way [31] . 

5 Summary and concluding remarks 

In this paper, we have discussed the basic notion of the nucleon-nucleon potential and its field-theoretical 
derivation from the equal-time Nambu-Bethe-Salpeter wave function in QCD. By construction, the non- 
local potential defined through the projection of the wave function to the interaction region (the inner 
region) correctly reproduces the asymptotic form of the wave function in the region beyond the range 
of the nuclear force (the outer region). Thus the observables such as the phase shifts and the binding 
energies can be calculated after extrapolating the potential to the infinite volume limit. Non-locality of the 
potential can be taken into account successively by making its velocity expansion, which introduces the 
velocity-dependent local potentials. The leading-order terms of such velocity expansion for the nucleon- 
nucleon interaction are the central and the tensor potentials. 

To show how this formulation works, some results in the quenched and (2+l)-flavor lattice QCD 
simulations are shown for relatively heavy pion masses, ~ 400, 500, 700 MeV. We found that the NN 
potential calculated on the lattice at low energy shows all the characteristic features expected from the 
empirical NN potentials obtained from the experimental NN phase shifts, namely the attractive well 
at long and medium distances and the repulsive core at short distance for the central potential. As for 
the tensor potential obtained from the coupled channel treatment of the ^Si-state and the '^Di-state, we 
found appreciable attraction at long and medium distances. 

As the quark mass decreases, the repulsive core and attractive well in the central potential, and the 
attractive well in the tensor potential tend to be enhanced. To make the deuteron bound state, however, 
it is necessary to go the lighter quark masses. We have also shown that the derivative expansion in terms 
of the local and energy-independent potentials works well at low energies for at least the quark masses 
studies above. 

There are a number of directions to be investigated on the basis of our approach. Among others, the 
most important direction is to carry out (2-|-l)-flavor simulations with a large volume (e.g. L = 6 fm) 
at physical quark mass (m^ = 135 MeV) to extract the realistic NN potentials. This will be indeed 
started soon as a first priority simulation at 10 PFlops national supercomputer "KEI" which will have 
full operation in 2012 at Advanced Institute for Computational Science (AICS) in Kobe, Japan [32] . 
Simulations of the three or more nucleons on the lattice are also a challenging problem to be studied in 
relation to the attractive binding of finite nuclei and to the repulsive effect in high density matter relevant 
to neutron stars. Study along this line has been recently started [33] l34]. 

If it turns out that the program described in this paper indeed works in lattice QCD with the physical 
quark mass, it would be a major step toward the understanding of atomic nuclei and neutron stars from 
the fundamental law of the strong interaction, the quantum chromodynamics. 

Note added: Recently, full QCD simulations of the hyperon potentials in the flavor SU(3) limit reported 
in [25] and discussed in Sec. |4|were extended to the lattice sizes L ~ 3 and 4 fm for the pseudo-scalar 
meson mass of 673-1015 MeV. By solving the Schrodinger eqaution with the flavor-singlet potential, a 
bound i7-dibaryon with the binding energy of 30-40 MeV was found [35] . Since the binding energy turns 
out to be insensitive to the quark masses, there may be a possibility of weakly bound or resonant H- 
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dibaryon even in the real world with lighter quark masses and with the flavor SU(3) breaking. To make a 
definite conclusion, however, the (2+l)-flavor lattice QCD simulations for i?-dibaryon with AA-A^S-EE 
coupled channel analysis is necessary. Such a direction is currently in progress [36]. See also a related 
recent work on iJ-dibaryon in a difi'erent approach [37] . 
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